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^ Abstract 

We explore the reconstruction of the gravitational lensing field of the cosmic microwave back- 
ground in real space showing that very little statistical information is lost when estimators of short 
range on the celestial sphere are used in place of the customary estimators in harmonic space, 
which are nonlocal and in principle require a simultaneous analysis of the entire sky without any 
cuts or excisions. Because virtually all the information relevant to lensing reconstruction lies on 
angular scales close to the resolution scale of the sky map, the gravitational lensing dilatation and 
shear fields (which unlike the deflection field or lensing potential are directly related to the observa- 
tions in a local manner) may be reconstructed by means of quadratic combinations involving only 
very closely separated pixels. Even though harmonic space provides a more natural context for 
understanding lensing reconstruction theoretically, the real space methods developed here have the 
virtue of being faster to implement and are likely to prove useful for analyzing realistic maps con- 

■^j- taining a galactic cut and possibly numerous small excisions to exclude point sources that cannot 

q be reliably subtracted. 
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I. INTRODUCTION 



Recently there has been much discussion of how the intervening mass from clustered 
matter between the last scattering surface (at z ~ 1100) and an observer at the present 
time distorts the CMB anisotropy by means of gravitational lensing [THH]. On the level of 
the two-point correlation function, this effect distorts the TT (temperature-temperature) 
correlation power spectrum j6] and mixes the EE and BB polarization power spectra as well 
as distorting them [THS]. Lensing also introduces non-Gaussianities that manifest themselves 
in the higher-point correlation functions [TQl E] • At the level of the three-point correlation 
function, to leading order there is no nonzero expectation value if we regard the lensing 
potential as a random field p21 [13]. However if we consider the CMB lensing potential as 
fixed, we find that expectation values of the form 

(T(W)> 4(0 (1) 

do not vanish, and this property may be exploited to recover or "reconstruct" the lensing 
field using estimators quadratic in T (or in E and B) [H]. 

Much effort has been devoted to developing an optimal reconstruction of the lensing 
potential in harmonic space, which implicitly assumes full sky coverage with no galactic cut, 
no bad pixels due to point sources that must be excised, and no nonuniform weighting to 
account for uneven sky coverage [ToTfiT] . For a reconstruction based on the temperature 
anisotropy alone, it has been shown how to construct the optimal quadratic estimator in 
this idealized context [18] , and the improvement that can be gained from using an even more 
optimal maximum likelihood estimator is marginal [19] , because the distortion due to lensing 
is small compared to the intrinsic cosmic variance and noise of the experiment (although 
this assumption is less valid at very large I for very clean maps where the lensing signal is 
dominant). For the exploitation of polarized anisotropics, the situation is somewhat more 
complicated. When the experimental noise is large compared to the B polarization mode 
generated by lensing, the situation is essentially the same as for a reconstruction using the 
temperature data [201 ET] . However at a higher sensitivity where the B signal is essentially 
entirely due to lensing, the quadratic estimator underperforms because the actual multipole 
moments rather than their averages should be used for the optimal weighting. In this case, 
the higher order corrections to the quadratic estimator present in the maximum likelihood 
estimator are no longer negligible [22J. 

In this paper we investigate a real space approach to lensing reconstruction. Under the 
ideal conditions often assumed and described above, this approach would naturally yield the 
same result as the conventional approach in harmonic space. Our interest however lies in 
considering slightly non-optimal estimators that have been modified to have a finite range 
so that cuts, excisions of pixels, and non-uniform coverage may be included in a simple 
and flexible way. We believe that such non-ideal but more robust local estimators defined 
in real space may prove superior for confronting the complications inherent in analyzing 
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real data [2~iH2"6] . Another advantage of the approach here is that the dilatation and pure 
shear provide separate and essentially independent lensing reconstructions which may be 
confronted with each other. This feature may prove useful as a way of diagnosing spurious 
signals, which are unlikely to affect the two reconstructions in the same way. Moreover the 
presence of two shear components enables one to estimate the noise of the reconstruction 
through the implied transverse displacement field, which is forbidden in weak lensing. 

Before proceeding to the details of this program, it is useful to consider the relation 
between the various descriptions of the lensing and the deformation of the anisotropies in 
real space. It is also useful to consider which angular scales contribute the most statistical 
weight to the lensing reconstruction. 

The lensing distortion of the CMB anisotropy on the surface of last scatter may be 
described in three ways: by a lensing potential $, by a deflection field £ = V$, or by the 
three components of the shear tensor 

V„V,<I>. (2) 

Even if we have simultaneous access to the entire sky, the descriptions $ and £ suffer from 
an ambiguity. $ cannot be distinguished from $ + (constant) and the vector field £ can be 
measured only up to a constant translation (or more properly a rotation in the presence of 
sky curvature). This is because if we know only the CMB power spectrum, a patch of sky and 
its translation necessarily have the same likelihood on account of isotropy of the underlying 
stochastic process. Consequently, the absolute translation due to lensing cannot be observed. 
By contrast, locally the shear and dilatation (which are gradients of the translation vector 
field) are completely well defined. This can easily be seen by considering the effect of a 
constant deformation described by a deformation matrix S relating the angular coordinates 
6, the actual coordinates on the celestial sphere of a point on the last scattering surface, to 
the coordinates 0' , the coordinates that the same point would have in the absence of lensing. 
1 We have 6' = SO where S = exp[— k]. (Note that we employ the flat sky approximation 
and assume that the deformation is small so that a linear treatment is adequate.) To linear 
order, the power spectrum is modified in the following way by this linear deformation, which 
preserves the homogeneity but not the isotropy of the underlying statistical process: 




C(\£\) -> C{£) 



, , d(HC(£)}) | 2 \ , (k + {£\-£1)+k x {2£ x £ 2 )\ d(ln[C(£)}) 



d(ln[£\) J V P ) d(hi[£]) 



(3) 



For the case of perfect scale invariance (i.e., a power law of the form C(£) oc £ 2 ) there is no 
change in the correlations due to the dilatation component of S, and similarly for a perfect 

1 In the sequel we shall, unless otherwise indicated, employ the flat sky approximation where the vector 6 
represents a point on the flattened celestial sphere and t represents a wavevector. At times summations 
over (l,m) shall also be used. 
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white noise spectrum (i.e., a power law of the form C{£) oc £°) there is no sensitivity to the 
pure shear components k+ and k x in the anisotropic (m = ±2) correlations. 

To the extent that the shear-dilatation components are slowly varying, we may construct 
estimators of k , k+, and follows 

k = No" 1 J d 2 9 J d 2 9' [T(0) T(0')-(T(0) T(6')) K=0 ] 

x J eMii - {e -° )] m) + v torn ) 

= No" 1 A' 1 [ d 2 9 [ d 2 9' [T(0) T{0') - (T(0) T(6')) K=0 \ K o (0 - 0') 

J A J A 

= ^J^O [T{0) {F. T){0) - (constant)] (4) 
and similarly 

= AT +jX - 1 f d 2 9 [ d 2 9' T{0) T{0') 

J A J A 

x fA ranW , C(£) d(ln[C (£)}) (cos(2#) 

X ](2nf em%t '^ U )l [C{£) + N{£)f d(\n[£}) ^sin(2tf) 

where the normalization factors are given by 

N = A f dH jC(£)] 2 ( d(\n[C(£)]) V 
^7 {2nY\C{£) + N{£)Y \ d{\n[£\) + J ' 
A f d 2 £ [C{£)f f d{\n[C{£)}) V 



N+ Nx 2 J {2^\C{£) + N{£)Y V d{\n[£\) ) (6) 

and A is the area of the domain. Here N(£) is the noise of the experiment being considered 
and serves as a cut-off at large £, above which there is very little exploitable information 
because of the low signal-to-noise of the CMB maps. 

The above expressions assume a spatially flat, two-dimensional domain formally of large 
but finite area subject to a spatially uniform linear deformation. We may consider a toroidal 
domain in the limit that the periods of the torus become arbitrarily large (though the 
toroidal domain is not essential for the application of the real space estimator presented in 
this paper). Before proceeding it is useful to consider the relation of this simplified problem 
to the real problem of reconstructing the lensing field starting from a CMB map of finite 
resolution. Because of the smallness of the lensing distortion of the CMB anisotropy, the 
idealized situation considered above is less different from the actual situation, where one is 
dealing with a curved sky and a varying dilatation and shear fields, than one might at first 
sight suppose. 
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FIG. 1: Lensing power spectrum. The lensing field power spectrum is shown represented in 
several manners. The three panels (from left to right) illustrate the lensing field expressed as po- 
tential, a deflection field, and dilatation/shear field, respectively. Plotted are [£{£+\)Cxx I '(2tt)] 1 ^ 2 
where XX = $$, kk. Here cf* = l{l + 1) Cf * and C?" = £ 2 (£ + l) 2 Cf */4. Of the three, 
the last power spectrum is more directly related to the observed distortion. 

The rightmost panel of Fig. [T] shows the shear-dilatation power spectrum as a function of 
multipole number £, and we observe that for £ < 100, the distortion is always less than about 
1.5%. This implies that in order to attain an (S/N) of approximately unity it is necessary to 
consider a region containing at least O(10 3 ) resolution elements, where a resolution element 
is a pixel of the map of sufficient size so that the noise and angular resolution of the survey 
give S/N « 1. Consequently, there is little point to trying to reconstruct the lensing field 
over a region not having at least 30 resolution elements on a side. If the distortion from 
lensing were greater, the situation would be different. 

For the ideal linear estimator 



K 0, ideal 



£ 

£,m 

A 

(47T) 



CI 



2(C t + N t y 



d(\n[C}) 
_ d(ln[£}) 

C(£) 2 



(C(£) + N(£)) 2 



d(\n[£}) 



1 2 



+ 2 



(7) 



where A is the area of the sky patch considered. 

We now consider other unbiased estimators of k,q and n + where a slight increase in the 
variance is compensated for by other desirable properties. We are presently interested in 
estimators for which the real space filtering kernel is short range. To this end it is useful to 
define the inner product on the space of weight vectors w = {we}: 



A 



2(C e + N t 



\2 W £,rn W £,m 



; w A (£)w B (£) 



(2ir) 2 2[C(£) + N(£)] 



(8) 
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where we give both the spherical and flat sky continuum forms. If we set Wid ea i{£) 



SC{e) th)K0=1 = C + 2], then in terms of the above inner product 



and 



\'> nl< ail SCobs) 

K 0,ideal — 7 T l y J 

{Wideal, W ideal ) 



(5) = ($Cth,K =u8C th;Ko=1 ) . (10) 

\ / Km, ideal 



kq, ideal 

Given an arbitrary weight vector w, using the above inner product we may define the fol- 
lowing unbiased estimator of k 

A (w,5C obs ) 

«oM = t r 11 

{W, W idea i) 

provided that the denominator does not vanish, and its variance is given by 
(w,w) I (wideai,w) 2 , so that the increase in variance with respect to the ideal estimator 
is given by following geometric expression for the secant squared 

V&y(k {w)) (W,W) (W ideal, Wideal) 2/ \ n \ 

sec [x). (12) 



V&r(k (w ideal )) (W, Wideal) 
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For k + and n x analogous formulae may be derived straightforwardly. 

We consider how the information (or S 2 /N 2 ) contained in the ideal estimator is dis- 
tributed over the various multipoles. In Fig. [2] [panels (c) and (d)] we plot the quantity 





If C(lf 


~d(hx[C(i)]) | o" 


2 




lt (C(£)+N(£)) 2 






Jo 00 ' 


di C{£)2 


\d(HC(i)}) o 


2 


(C{£)+N(£)) 2 







F < = , 7 , , ' 2 (13) 



and F>(£) = 1 - F<(£) where 

iV £ = iV exp [+£ 2 e 2 b eam] = N exp [+f/£t am ] (14) 

and ^beam = (810)(10V^6eanT)- The corresponding quantity is also shown for the shear. 

In their present state, the ideal minimum variance kernels for the estimators k , k + , and 
k x have their support sharply peaked at small separations, but nevertheless there is still 
some small support for large separation. We now investigate how much information is lost if 
the support at large separation is completely cut away. The quantitative way to characterize 
this loss is to ask by what factor the variance of the estimator is increased relative to the 
optimal estimator after our pruned estimator has been renormalized to render it unbiased. 

In real space the minimum variance full-sky estimator kernel for kq is given by 

POO 

KidealiO) = \ tdi K idea l{l) 
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FIG. 2: Character of the signal. The top row shows as a function of multipole number £ the 
temperature power spectrum and local spectral index, defined as dlnCi/dlni, for the standard 
cosmology (WMAP best-fit model). On the bottom row, the left panel shows the normalized 
cumulative x 2 as a function of £ integrated both from the left and from the right using the sensitivity 
and resolution parameters for the PLANCK experiment (where the 100, 143 and 217 GHz channels 
have been combined in quadrature) [23J. We observe that the central 80% of the information is 
concentrated in the range £ = 800-1600. Smaller £ contribute almost no information because there 
are comparatively very few independent multipoles, and moreover the angular spectrum in the sky 
is very nearly scale invariant. At much larger £ instrument noise and beam smearing wash out 
the usable signal. In the intermediate range a structure of plateaus connected by steep rises can 
be observed. This structure is a direct result of the Doppler oscillations. Around the crests and 
troughs the spectrum is almost scale invariant and hence does not contain any information for 
determining the dilatation. The right panel shows the corresponding plot for the shear, where the 
plateaus are less pronounced. 



idi J (£9)- 1 ( ! '' ! 



d[\n(£)] 



+ 2 



(15) 



N ideal [C{£) + N{£)Y 
This kernel may be inverted using the following inverse Bessel transform: 

1 f°° 

K ideal (£) = — 6d6 J Q (£6) K ideal {9). (16) 
2vr Jo 

We limit the support of the kernel by requiring that K{9) be nonzero only for 9 < 9 max where 
9 max is varied. This is accomplished numerically by expressing Jq(£9) as a linear combination 
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FIG. 3: Performance of estimator with truncated angular support. We indicate how 
limiting the angular support of the filter in our estimator increases the noise. The panels on the 
top row refer to the dilatation filter, while the panels on the bottom row refer to the shear filter. 
Panel (a) indicates how the estimator variance (with the estimator normalized to be unbiased) 
increases as the angular support (disk radius in degrees) is reduced. Panels (b) and (c) indicate 
the profiles of the optimal truncated estimators in both angular space and harmonic space. 



of cubic spline basis functions spanning the interval 9 e [0, 9 ma x\ and optimizing for the shape 
that minimizes the variance calculated according to eqn. (12). Analogous expressions may 
be obtained for the shear by replacing J with J 2 . Fig. [3] shows the variance ratio as a 
function of 6 max for the dilatation and shear estimators. We observe that the increase in 
variance at small separations is more modest for the shear estimator. 



II. RESULTS: DILATATION AND SHEAR RECONSTRUCTION 

In the previous section we showed how by means of a linear filter T K0 applied to a 
temperature map T(0), we may obtain a reconstructed dilatation field through the product 
map 

K 0irec (0)=T(0)(T Ko T)(0)-c (17) 

where c is a constant offset. We presented a theoretical derivation of the optimal shape for 
such a filter. The best shape for the filter depends both on the cosmological model, and 
more importantly on the details of the experiment, because most of the statistical weight is 
situated on the smallest angular scales, near the resolution limit of the experiment. If the 
primordial anisotropies dominated to arbitrarily small scales and experiments of unlimited 
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sensitivity and angular resolution were possible, the lensing field could be reconstructed 
as accurately as desired simply by increasing the angular resolution and sensitivity. For a 
given experiment, the desirable filter shape may be intuitively understood as follows. We 
assume weak lensing (i.e., \k\ <C 1). We want to measure the change in shape of the power 
spectrum (or rather of the "local" power spectrum within a certain finite size patch for the 
case of interest of a non-uniform dilatation). Working in the flat sky approximation, which 
is appropriate because the length scales most sensitive to the changes in overall scale are 
small (i.e., £ ^> 100), we define C(£) = £ 2 C(£). Under a small dilatation, which transforms 
the temperature map as follows 

T(0) T{0') = T(exp[-K] fl), (18) 

the power spectrum transforms as 

C(£) ->c(exp[+K] (19) 

so for positive dilatations the power spectrum in the form C{€) squeezes to the left in such 
a way that the amplitudes of the various features remain unchanged. If the dilatation is 
small, the change in the power spectrum is proportional to k C'(£). To make a good filter for 
implementing the above scheme giving the reconstruction with the least noise, we wish to 
block modes of wavenumber where C(£) is approximately flat, to allow modes of wavenumber 
where C(£) is rising to pass with a positive phase factor, and to allow modes of wavenumber 
where C(£) is falling to pass with a negative phase factor. Moreover the filter should combine 
the different wavenumbers according to inverse variance weighting in order to minimize the 
noise of the reconstructed dilatation field. 

The quality of the inverse variance reconstruction may be characterized quantitatively in 
terms of the variance of the reconstructed field. We define a quality factor Q so that 

(8Kf » ? (20) 

where A is the area of the region over which the reconstructed field has been averaged and 
(5k) is the fluctuation in k averaged over that same area. Here Q is equal to a\ Q / 'A where 
cr? o is given in eqn. (|7j). 

Implicit is the assumption that the noise in the reconstructed dilatation field lacks long 
range correlations — that is, it can be treated as white noise beyond a certain angular coher- 
ence scale. This is indeed the case because the filter J-" Ko blocks small wave numbers. We 
assume that the area A is sufficiently large so that the white noise regime has been reached. 

For a non-uniform dilatation field having a small wavenumber, the previous analysis 
carries over with very small corrections. Using the same filter as for the uniform case, one 
expects that 

(ftO,recon(0)) = K O,exarf(0); (21) 
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FIG. 4: Estimator form factors at large lensing wave numbers. We plot as a function 
of multipole number the form factors for the dilatation, longitudinal shear, and cross shear re- 
constructions using the filter described in the text. Noise for the Planck instrument combining 
the 100+143+217 GHz channels is assumed using bluebook values and combining channels in the 
naive way |23j . A deformation field whose wavenumber constitutes the horizontal axis was applied 
to unlensed maps. For the dilatation and the "plus" shear, a longitudinal displacement field was 
used. For the "cross" shear a transverse deformation was applied. As a result of symmetry consid- 
erations, the responses of the dilatation and plus shear to a transverse deformation have vanishing 
expectation values (i.e., they are pure noise) and hence are not shown. The same holds for the 
cross shear response to a longitudinal deformation. 



in other words, that there is negligible bias. However for larger wavenumbers a form factor 
appears due to the fact that our estimator probes the power spectrum over a window of a 
finite width in angular space, and consequently smoothes the underlying exact dilatation 
field thus reducing its amplitude. Quantitatively, if the dilatation field has the form 



k cos(£ ■ + 0), 
the expectation value of the reconstructed field takes the form 



(22) 



F(£)k cos(£-6 + 



(23) 
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FIG. 5: Alternative filter derivation. We generate a dilatation reconstruction filter by an 
alternate procedure whereby Gaussian niters are placed at the locations in multipole space where 
the magnitude of the derivative of I 2 C(l) is greatest. The widths of the filters are chosen to 
correspond roughly to the width of the rise and fall of this quantity. 

and as before the noise field will have the same statistics. F(£) is known as the form factor 
and satisfies -F(O) = 1. F(£) falls off for increasing £, eventually going to zero. The form 
factors are indicated in Fig. |4} The shear plus estimator experiences less smoothing than 
the dilatation estimator because the form of the ^-space filter has less cancellations as can 
be seen in Fig. |3j 

In Fig. [5] we present an alternative, more heuristic derivation of the dilatation field re- 
construction filter which has the virtue of rendering manifest the contribution of various 
features in the power spectrum to the reconstruction. This alternative filter is a combina- 
tion of Gaussian filters placed at the locations in multipole space where the magnitude of 
the derivative of £ 2 C(£) is greatest, and whose widths are chosen to correspond roughly to 
the width of the rise and fall of this quantity. In Table [i] the (S/N) 2 for each of these filters 
are shown and we indicate their relative statistical weights as well as the (S/N) 2 when they 
are combined using inverse variance weighting. Somewhat surprisingly, this naive derivation 
yields a filter almost as efficient as the optimal one. 

We now repeat the same analysis for a constant pure shear deformation, for concreteness 
with k + > 0, k x = 0, there is a stretching along the x-direction with the same amount of 
compression along the y-direction, so that the temperature map is deformed according to 

T(0„ e y ) T'(6' x) e' y ) = T(e-^6 X , e + ^6 y ), (24) 
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3.13 ±0.26 
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1.56 ±0.13 
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40 


3.02 ±0.25 


770 


40 


0.13 ±0.01 


900 


50 


61.51 ± 5.02 


1070 


40 


1.44 ±0.12 


1200 


50 


30.66 ±2.50 


1370 


50 


23.24 ± 1.90 


1550 


40 


42.53 ±3.47 


Combined 


167.22 ± 13.65 


Ideal 


140.50 ± 11.47 



TABLE I: Performance of individual feature niters. The first nine entries correspond to 
Gaussian filters whose centers and widths are indicated and have been placed where the features of 
the power spectrum change the most in response to a small dilatation. The last column expresses 
the sensitivity of each filter. The next-to-last row shows (S/N) 2 resulting when the filters are 
combined using inverse covariance weighting. This is compared to the ideal filter. The errors result 
from Monte Carlo noise. 




FIG. 6: Distortion of power spectrum due to shear. In the two-dimensional plot in harmonic 
space the power spectrum is shown after a pure shear transformation has been applied turning the 
sequences of acoustic peaks into concentric ellipses rather than concentric circles. In the real space 
images (middle, unlensed and right, sheared) the transformation is hardly apparent to the unaided 
eye. 

and for the harmonic coefficients one has the following transformation law 

T(£ x ,£ y ) -> T'(£' x ,£' y ) = T(e +K +£ x ,e- K +£ y ). (25) 
For any power spectrum having a shape like that of the CMB temperature, where at all 



12 



wavenumbers the spectrum is "red" compared to a white noise spectrum, a pure shear 
deformation causes a loss of power for wavevectors oriented close to the stretched principal 
axis and a corresponding increase of power for wavevectors close to the principal axis of 
compression. Along the diagonal direction the power is unaffected. The estimator developed 
in the first section of this paper exploits this anisotropy to reconstruct two components of 
the shear field. The effect is illustrated in Fig. [6] In panel (a) the acoustic oscillations, which 
would in the absence of lensing be visible as a series of concentric circular rings, have as a 
result of a pure shear deformation been deformed into concentric ellipses. Panels (b) and 
(c) show simulated temperature maps before and after a constant shear deformation, here 
exaggerated in magnitude (with kq = 0.1) for clarity. One observes that the annuli become 
elliptical. Unlike for the dilatation field, where a scale-invariant spectral index implies that 
statistically the pattern does not change, for pure shear an isotropic scale invariant random 
field is deformed into an anisotropic scale invariant random field, where heuristically one 
might say that the series of elliptical (rather than on the average circular) diffuse motifs 
have been put down in a scale-free manner. In Fig. [7] we show the extent of the region of 
linear response of the dilatation and shear estimators. Ignoring higher order terms gives rise 
to a bias which can be corrected [HI 127] . 

III. COMBINING THE k (9), k+(6) AND k x (0) RECONSTRUCTIONS 

In the previous two sections we have shown how by means of three filters J- Ko , J- K+ , and 
J-" Kx we may obtain a noisy reconstruction of the fields Ko(0), k + (6), and k x (0) in real space. 
We also saw how the reconstruction was essentially finite-range in real space and could be 
made of compact support with very little loss of information relative to the optimal nonlocal, 
full-sky reconstruction. The reconstructed fields are not independent but related according 
to the following relations 



which are local in harmonic space. These consistency relations, expressed in terms of the 
lensing potential $, can be used to harmonize the reconstruction. Even though there are 
three reconstructed fields, the cross shear, which would result from a transverse deformation 
that cannot be produced by weak lensing, is pure noise. The noises present in the longitudinal 
shear and the dilatation reconstruction are uncorrelated because one involves the m = 
sector and the other the m = ±2 sector. Harmonization by inverse variance weighting 
reduces the noise. If we consider the noise in the reconstruction field to lowest order in 
the perturbation expansion (as 4-point functions of the underlying Gaussian field) ignoring 
higher-order terms, then there are no correlations between reconstructed fields components 



K + (£) 
K x {£) 



-(£l + f y M£), 

l(£l-e y M£), 

U y ®{£), 



(26) 
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FIG. 7: Estimator non-linearity. We plot the recovered root-mean-square distortion versus 
the input root-mean-square distortion in order to characterize the range of linear response for 
our estimator. The input is a long-wavelength longitudinal deformation, so that the two-dimension 
shear and dilatation are exactly half the one-dimensional dilatation. Exactly the same non-linearity 
will plague the quadratic estimator as well, since the two have been demonstrated to be equivalent 



at low wavenumber (as discussed in section IV) 



of differing wavenumber. Consequently, in order to obtain the lowest noise estimator of $, 
it suffices to consider each £ sector separately and to use the covariance matrix of the noise 
in Kq(£), and k x {£) as a basis for the optimal harmonization. 



IV. COMPARISON WITH THE QUADRATIC ESTIMATOR 

In this section we show how the dilatation and shear estimators developed in this paper 
are related to the linearly optimal quadratic estimator for the lensing field [T4"l ITS] . We 
demonstrate that in the low-£ limit the quadratic estimator becomes a linear combination 
of our dilatation and longitudinal shear estimators. In fact in this limit the linearly optimal 
quadratic estimator is identical to the inverse variance weighted linear combination of the 
dilatation and longitudinal shear estimators, whose noises are uncorrelated. At higher £, 



14 



the approximation used to derive this relation progressively breaks down. To quantify this 
degradation, we plot in Fig. [8] the increase in variance of the combined dilatation plus 
longitudinal shear estimators. For the reasons already presented concerning the minute 
statistical weight contributed by small multipoles, we find it adequate to work in the flat 
sky approximation. 

The relation in real space (valid to linear order) 



ST(0) = V$(0) • VT(0) 
in harmonic space translates into 



ST(£) 



d 



2 



(27T) 



'-£') ■ {£-£') T{£-£') 



where we define 

(T(£)T(£')} = (2n) 2 5\£ + £') C(£). 
Since T s k y = T + 5T, it follows that to leading (linear) order 
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is the minimum variance unbiased estimator for the Fourier coefficient $(£') where 
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We may approximate the quantity in the square brackets in eqn. (30) to linear order to 
obtain 
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where G is the angle between £ and £' . The expression is remarkably similar to a linear 
combination of the expressions appearing in the dilatation and the pure shear reconstruction 
of the previous section. 



Approximating the numerator of (32) as (34) and the denominator of (32) as [C(£) + 
N(£)] 2 , we may write the three estimators in the following unified manner: 
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and Nq = Ne> + Ns L because the cross term cancels under the integral. 

The quantities Nd, Ns L , Nq express the amount of information provided by the respective 
estimators, and we have shown that in the limit where £' is small 
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The exact expressions in Fourier space for the dilatation and shear statistics at not 
necessarily small wavenumber £' are 
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Multipole number (I) 



FIG. 8: Comparison with quadratic estimator We plot the noises of the various estimators 
compared to the expected signal (heavy red curve). The quadratic estimator is indicated in thick 
black. The dilatation and shear estimators are shown in dashed red and green, respectively, and 
when combined nominally give the dashed magenta curve, but when the imperfect overlap with 
the expected signal is taken into account, yield the solid magenta curve. The blue curve would be 
indicative of the actual noise in the recovered maps, but if the imperfect overlap where corrected 
to remove the bias at high £ the heavy magenta curve would result. For comparison we show the 
predicted lensing signal (as computed by CAMB for the WMAP best fit model) as the heavy red 
curve. 

where 
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To express compactly the increase 
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zero, it is useful to define the inner product 
(A,B 

where 
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and 5 is defined analogously. 

In this notation the various estimators have the following shape functions: 
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It follows that the increase in the variance is given by the following expression for the secant 
squared 

'Q,Q)(D.D 
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(43) 



In Fig. [8] we indicate the performance of the several estimators at finite £'. We plot 
(i.e. the expectation value of the square modulus of the harmonic coefficient on 
the sphere for the reconstruction noise) for the estimators Q, D, and Sl, in black, dashed red 
and dashed green, respectively. The dashed magenta curve indicates the noise obtained by 
combining D and Sl using inverse variance weighting. The solid magenta curve indicates the 
correction when the lack of overlap is taken into account to render the estimator unbiased 
at large wavenumber. We assume an instrument noise combining the Planck 100+143+217 
GHz channels according to the specifications given in the bluebook [23] ■ At low wavenumber, 
Q exhibits a flat (i.e., white noise) spectrum, which subsequently divergences at large £'. At 
very low £' we observe that the noise from Q is the same as the noise from D and Sl com- 
bined in quadrature, as shown theoretically in the text, but at higher £' the approximations 
used break down and a noise excess is observed. We observe that for £ <> 70, the difference 
in performance between the estimator developed here and the linearly optimal quadratic 
estimator is minimal. At higher wavenumbers beyond £ > 100, however, the variance in- 
creases rapidly due to lack of overlap with the ideal kernel. There is a priori no reason why 
a real space approach could not be extended to higher wavenumbers for the lensing field. 
However, in the present paper we do not explore how this would work. 
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V. CONCLUDING REMARKS 



We have demonstrated how to reconstruct in real space using a filter of compact support 
the weak gravitational lensing field, here represented as three fields, a dilatation field kq(0) 
and the two components of the pure shear distortion field k+(0) and k x (0). The three 
fields are related by a set of nonlocal consistency conditions, which may subsequently be 
exploited to reduce the noise of the reconstruction. Except for an integration constant 
and two translational and one rotational zero modes, the weak lensing may alternatively 
and equivalently be described by either (1) a gravitational lensing potential (2) a 

displacement field $(0), (2) a displacement field £(0), or (3) the dilatation field Kq(8) 

and the two components of the pure shear distortion field K+(6) and n x {0)- In this paper we 
argue that for the purpose of reconstruction the representation (3) is advantageous because 
this is the representation for which the lensing field bears a local relation to the real space 
CMB maps. This locality comes at a price because the three components are not independent 
and subject to nonlocal consistency conditions, which may be exploited to improve the 
reconstruction. Locality allows different regions of the sky to be analyzed independently in 
a natural way, quite unlike the quadratic reconstruction in harmonic space, where the entire 
sky must be analyzed simultaneously. This approach and variations thereof hold promise 
for dealing with partial sky coverage and excised point sources. For the filters developed in 
this paper there is very little loss of information for the lensing field at low wavenumbers. 
However at larger wavenumbers the lensing signal is attenuated according to a wavenumber 
dependent form factor, which can be deconvolved by applying a correction filter. 
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